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DERIVATION  OF  THE  BOLTZMANN  EQUATION  AS 


A  TEST  CASE  OF  KINETIC -THEORETICAL  SCHEMES  ♦ 

by 

Toyoki  Koga** 

Polytechnic  Institute  of  Brooklyn,  Graduate  Center 
Farmingdale,  New  Y  ork 

SUMMARY 

The  derivation  of  the  Boltzmann  equation  from  the 
Liouville  equation  can  be  a  means  of  testing  the  feasibility 
of  a  kinetic -theroeticad  scheme,  since  the  Boltzmann 
equation  is  known  as  empirically  plausible.  From  this 
viewpoint  •  we  investigate  kinetic -theoretical  schemes  deriving 
the  Boltzmann  equation.  We  may  derive  the  Boltzmann  equation 
consistently  by  considering  the  Liouville  equation  to  govern  the 
distribution  of  a  system  in  the  phase  space  which  is  not  symmetric 
in  the  microscopic  sense.  On  the  other  hand,  we  cannot  derive  the 
equation  consistently  on  the  assumption  of  symmetry  of  the 
distribution;  the  Boltzmann  collision  integral  is  not  compatible  with 
the  symmetry  assumption. 
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Agency  under  Order  No.  5^9  through  the  Office  of  Naval  Research. 

Visiting  Professor;  currently  on  professional  staff  of  TRW  Systems, 
Inc. ,  Redondo  Beach,  California. 
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I.  INTRODUCTION 


By  examining  usual  procedures  of  the  derivation  of 
the  Boltzmann  equation  from  the  Liouville  equation,  we 
find  inconsistencies  in  the  procedures.  The  main  purposes 
of  this  paper  are:  1)  to  point  out  a  common  defect  in  the 
functions  of  modern  kinetic -theoretical  schemes  by 
analyzing  these  inconsistencies  and  1)  to  att'^mpt  eliminating 
the  inconsistencies;  the  attempt  is  the  one  to  rationalize  kinetic 
theory  itself. 

The  Boltzmann  equation  is  a  product  of  intuitive  but 
ingeiious  thoughts  of  Maxwell  and  Boltzmann.  *  In  spite  of 
the  intuitiveness  involved  in  the  original  derivation,  the 
feasibility  of  the  equation  has  been  tested  and  confirmed 
in  many  ways  ^  ;  the  equation  is  to  be  appreciated  as  a  basic, 
if  empirical,  law  by  itself.  Therefore,  as  a  means  of  proving 
the  feasibility  of  the  Boltzmann  equation  itself,  it  is  of  no 
significance  to  show  a  radical  derivation  of  the  equation  from 
the  Liouville  equation.  On  the  contrary,  it  is  more  likely  to 
be  a  proof  of  the  feasibility  of  a  kinetic -theoretical  scheme  that 
the  scheme  is  able  to  derive  consistently  the  Boltzm2Uin  equation 
from  the  Liouville  equation.  In  other  words,  the  derivation  of 


^  It  is  well-known  that  the  equation  was  proposed  by  Boltzmann 
in  187^.  We  may  see,  however,  that  Maxwell' s  equations  of  trans¬ 
fer  proposed  in  1866  would  have  stimulated  Boltzmann  to  conceive 
the  equation.  See  The  Scientific  Papers  of  James  Clark  Maxwell 
(Dover, N.Y.),  Vol.  H,  p.  ib. 

^  S.  Chapman  and  T.  G.  Cowling,  The  Mathematical  Theory  of  Non  - 
Uniform  Gases  (Cambridge  University  Press,  1^61),  second  edition. 
For  example,  by  assuming  a  proper  law  of  interaction  force  between 
molecules  constituting  a  gas,  we  may  oalculate  the  viscosity  co¬ 
efficient  and  heat  conduction  of  the  gas  according  to  the  Boltzmann 
equation. 
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the  Boltzmann  equation  under  certain  conditions  can  be  a 
trial  of  the  feasibility  of  kinetic-theoretical  scheme.  i  Of 
course,  there  may  be  a  benefit  that  such  a  derivation  shows 
us  more  explicity  conditions  which  limit  the  validity  of  the 
Boltzmann  equation. ) 

There  have  been  attempts  made  by  several  authors  to 
derive  the  Boltzmann  equation  from  the  Liouville  equation. 

In  spite  of  their  radical  intention,  however,  the  attempts  do 
not  seem  suoL^ssful:  For  conceiving  properly  the  Boltzmann 
equation,  it  is  necessary  to  make  the  concept  of  molar  disorder 
and  that  of  molecular  order  as  compatible.  The  compatibility 
of  these  two  conflicting  concepts  of  states  of  particles  in  the 
Boltzmann  equation  was  made  possible  ingeniously  by  Boltzmann 
and  Maxwell  by  means  of  the  Stosszahlansatz.  On  the  other 
hand,  by  a  purely  kinetic -theoretical  approach,  the  compatibility 
tends  to  be  made  possible  only  on  the  destruction  of  consistency 
of  the  analyses.  The  common  failure  of  known  kinetic -theoretical 
schemes  deriving  the  Boltzmann  equation  is  found  at  this  point. 
Then,  changing  our  viewpoint,  the  very  failure  suggests  to  us  that 
there  exists  a  basic  and  common  defect  in  the  structures  of  these 
kinetic -theoretical  schemes.  That  is  the  assumption  of  the 

See  section  4. 


microscopic  symmetry  of  the  distribution  of  a  system  in  the  phase 
space*.  Based  on  this  as*.umption,  modern  kinetic -theoretical 
schemes  are  not  well-equipped  with  machinery  by  which  it  is 
possible  to  treat  molar  disorder  and  molecular  order  aj 
compatible. 

In  section  2,  basic  concepts  and  definitions  necessary 
for  the  dervation  are  introduced.  It  is  stressed  that  the 
distribution  of  a  system  in  the  phase  space  is  not  symmetric 
in  the  microscopic  sense.  In  section  3,  it  is  shown  that  coarse - 
g-.'aining  operations  with  respect  to  time  and  also  with  respect 
to  similar  particles  are  necessary  for  deriving  the  Boltzmann 
equation;  the  compatibility  between  molar  disorder  and  molecular 
order  is  assured  by  the  difference  between  the  macroscopic  and 
microscopic  time  scales.  In  comparison  with  the  derivation 
given  in  section  3,  we  identify  the  common  difficul^y  and  its 
causes  in  known  derivations  by  taking  the  Kirkwood-Ross 
derivation  and  the  Bogoliubov  derivation  for  examples  in 
section  4.  We  also  evaluate  the  derivation  proposed  by  Grad 
which  is  rather  anomalous. 

e 

The  distribution  is  not  symmetric  in  the  microscopic  senne. 

T.  Koga,  "Reply  to  Comments  of  P,  Schram,  Phys.  Fluids, 

6,  (1963);  "A  Kinetic  Theoretical  Investigation  of  a  Fully 

Ionized  Gas",  Part  1,  PIBAL  Report  No.  863,  January  1965. 


II.  PREPARATORY  DEFINITIONS  AND  TREATMENTS 


The  Liouville  equation  governing  the  evolution  of 
the  distribution  of  a  system  consisting  of  N  similar  particles 
(nriaterial  points)  is 


( 


D^'^^  =  0 


(2.  1) 


with 


=  V  =  -Jl.,  =  0 


U.2) 


( N) 

Here  D  '  is  the  distribution  function  of  the  system  in  the 

phase  space,  t  is  time;  m  is  the  mass  of  a  particle;  p^  is 

an  independent  variable  representing  a  point  in  the  momentum 

space  for  particle  i;  r^  is  an  independent  variable  representing 

a  point  in  the  position -vector  space  for  particle  i;  ^  is 

the  total  force  exerted  on  parlhcle  i;  .  is  the  force  exerted 

^  J 

on  particle  i  by  particle  j.  It  is  assumed  that  there  is  no 

external  force  exerted  on  the  system.  Since  we  are  investigating 

(N) 

a  single  system,  D  is  to  be  given  by 


D 


(N)  ^  N 
i=  1 


6(X.  -  X.  (t)) 


(2.3) 
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where  is  a  six-dimeneional  .ctor  standing  for  p.  and  r. 


X,  =  p..  r. 
1  1 


and  (t)  constit'ites  a  set  of  solutions  of  the  equations  of 

mot ’on  of  the  N  particles.  It  is  noted  that  is  not 

symmetric  with  respect  to  the  interchange  of  coordinates 

* 

between  any  two  particles  : 


x.«  •  •  •  I  X,, 
»  J 


Xj . X., 


Xj^)  (2.4) 


Distribution  function*  of  subsystems  are  defined  as  follows: 


(X.;t)  =  V  J  n  dX. 
'  '  iti  ^ 


(2.5) 


where 


dX.  =  do.  dp.  dp.  dr.  dr.  dr. 

J  *jx  J*  jy  J* 


and  V  is  the  volume  of  a  domain  of  the  configuration  space  in 
which  the  system  is  known  to  exist  throughout  the  entire  period 
of  time  of  our  investigation;  the  domain  may  be  chosen  arbitrarily 
as  long  as  it  is  sufficiently  large  that  the  system  or  its  part  will 
not  run  off  the  domain.  Once  it  is  chosen,  however,  the  domair. 


See  footnote  on  page  4. 
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should  not  be  changed  during  our  investigation.  Further  we  define 

(2.6) 


F  ;X  ;t)  =  n  dX,  .  if  i  ^  j 

•’  J  k  i,  j  ^ 


=  0  ifi  =  j  (2.6)' 

and  so  forth.  It  is  notod  tha.t  the  order  of  subscripts  attached  to 

(2) 

F'  '  is  the  same  as  the  order  of  t'  e  same  set  of  subscripts  at¬ 
tached  to  the  X' s  in  F^^^.  This  accordance  of  the  orders  means 
that  there  has  been  no  interchange  of  coordinates  among  particles. 
According  to  the  definition,  it  holds  that 


F<?/(X..  X.;  t)  =  F.  .<2>(X..  X,;t) 


(2.7) 


and 


■jl*'  {Xj.  X.;  t)  i  Py'^’(Xj,  Xji  t)  (2.7)' 

in  view  of  (2.  4). 

By  ii  tegrating  Eq.  (2.  1),  term  by  term,  with  respect  to 
the  coordinates  of  all  the  particles  except  for  particle  i,  we  obtain* 


P. 


'-It  +  ^  •  -3^  1 


(continued  on  next  page) 


Because  of  the  finiteness  of  V  and  of  the  energy  of  the  system, 
it  holds  that 


Therefore 


/4  °  "A;  “ 


(N) 


dr,  =  0 
k 


/ 


F,.  .  F  dX,  =  0 

)•  ^  J 


This  term  is  retained  in  equation  (2.  8)  for  the  convenience  of  the 
further  treatment  of  the  equation. 
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(2.  8) 

If  the  integration  is  ma<’.e  with  respect  to  all  the  particles  except 
for  particle  i  and  particle,  j.  we  obtain 


(2.  9) 


We  may  ir4crease  the  number  of  particles  of  exception  so  that 
the  equations  of  evolu'J  -n  of  and  others  are  derived  successively. 
By  making  averages  with  respect  to  particles,  we  define 


U) 


OJ 


I 

i 


X  +  aX 

F.‘  \k..  X..  t)  dX. 


F^^hx 


x^-.t) 


A  X  AX, 
a  D 


(2.  12) 


r 

i  j 


X  +  aX. 

I  F\/(X,  X.;t;  dX.dX. 
,j  'J  "  J  1  J 
b 


U.  13) 


Functions  of  X.  and/or  X^  are  discontinuor.j  functions  as  shown 
by  (2.  3).  but  functions  of  X^  and/or  X^  may  he  continuous  as 
the  number  density  of  particles  increases.  In  the  above,  it  is 
worthwhile  to  recall  that  X^,  Xy  etc.  are  independent  variables. 

Defined  by 


f(c,  r;t)drdp/m3  =N  F^^)(p.  r;t)drdp 


(2.  14) 


c  =  p/m  , 


dr=dr  dr  dr  , 
X  y  z 


etc. 


f  is  the  number  density  of  the  particles  specified  with  velocity 
c  and  position  vector  r  at  time  t. 

It  ii  often  aeaumed  that  f  defined  by  (2.|4)  ie  the  very 
function -yhich  evolve,  according  to  the  BoU.mann  equation.  Wt 
note,ho„ever,  that  there  i.  a  difference  between  f  and  the  function 
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to  be  governed  bv  the  Boltzmann  equation  as  is  e:  -  lained  in 
the  following.  The  Boltzmann  equation  is  postulf  ted  by 

+  c.^)£'{c.r;t) 

=  "  [  f  (c\ ,  r;  t)  f  ic' .  r;t)  -  f  (c,  r;t)] 

|c, -  c|  dy  hdhdc^  (2.15) 

According  to  the  definition  of  the  Boltzmann  equation,  c,  c, ,  c*  and 
c'^  must  be  velocities  of  particles  only  in  free  flight;  therefore 
f  should  be  defined  so  that  f  gives  the  number  density  of  such 
particles  only.  Angle  to  length  h  are  variables  specifying 
the  geometrical  relation  of  a  collision  between  Iwo  particles  which 
begins  with  initial  (asymptotic)  velocities  c  and  c  ^  of  the 
two  particles,  as  are  illustrated  in  Fig.  2.  1;  c'  and  c’ ^  are  the 
final  (asymptotic)  velocities  of  the  two  particles  after  the  collision. 
In  Eq.  (2.  15),  c,  Cj,  and  h  are  independent  variables;  in  accordance 
with  the  equations  of  motion  of  two  particles,  a  collision  may 
begin  either  with  c  and  c^  or  with  c'  and  c',  and  ends, respectively 
either  with  c'  and  c'^  .  By  the  Boltzmann  equation,  the  evolution 
of  f*  is  due  to  the  effects  of  many  collisions  accounted  for  in  such  a 
manner  that  the  effect  of  each  collision  is  the  integral  of  the 
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relevant  differential  effect  over  the  entii  e  phase  of  the  collision. 

On  the  other  hand,  the  evolution  of  according  to  Eq.  (1.  8) 

ia  due  to  the  effects  of  collisions  accounted  for  at  given 
which  may  be  a  state  of  particle  i  during  a  collision  as  well 
as  during  a  free  flight.  See  the  schematic  comparison  given  in 
Fig.  2.  2.  Therefore  f  defined  by  (2.14),  which  counts  in  net  only 
the  number  density  of  pax  'ides  in  free  flight  but  also  that  of 
particles  during  their  collisions, is  not  the  very  distribution 
function  to  be  governed  by  the  Boltzmann  equation.  However, 
if  the  portion  of  a  unit  time  during  which  a  particle  is  in  interaction 
with  others  is  much  smaller  than  the  other  portion  during  which 
the  particle  is  in  free  flight,  we  may  eliminate  small  and  sudden 
fluctuations  in  f  due  to  r>articles  in  collisions,  by  making  anaverage 
of  f  over  a  time  period  from  t  to  t  +  t  . 


<'>o  =  T-  J 


f(c,  r(s);  t  s)ds 


(2.  16) 


where 


=  c.  r(0)  =  r 


If  T  and  r(T)-r  are  negligibly  small  in  the  macroscopic  sense 


. . 


and  yet  large  in  the  microscopic  sense.  <f>Q  defined  above  is 
equivalent  to  f  in  the  macroscopic  sense.  In  the  microscopic 
sense,  <f>Q»  unlike  f  itself,  is  a  smooth  function  of  t  presenting 
only  the  density  of  pr.rticles  in  free  flight.  Therefore  it  is 
more  plausible  to  take  <f>Q.  instead  of  f.  for  the  distribution 
function  f  to  be  governed  by  the  Boltzmann  equation  (2.  15)  ^ 

f  =  <f  Q  (2.  17) 
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m.  DERIVATION  OF  THE  BOLTZMANN  EQUATION 


First,  we  assume  that  the  significant  interactions  among 
particles  are  always  binary.  In  other  words,  if  is 

significant,  then  j;  and  /or  are  negligible.  According 

to  this  assumption,  the  integral  term  in  Eq.  (2.  9)  vanishes: 


'rvt  m  3r^  m  ^rJ 


+  +jr.  )  fA^^X.,  X.;t) 


'ij  ^Pi 


=  0 

(3.1) 

Secondly, we  assume  that  the  scale  Tj  of  gas  dynamical  processes 
in  tiie  system  is  much  longer  than  the  time  scale  t,  of  inter - 
particle  interaction: 


Tj  >■>  X, 

During  the  time  period  from  t  to  t  +s,  we  consider  fictitious 
trajectories  of  particle  '  which  satisfy 


dr.(s)  p. 

— 4 i —  ,  r.(s)  =  r.  +/jr  .(s) 
ds  m  i'  '  i  **  i'  ' 


(Trajectories  0) 


(3.2) 
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with  initial  conditions 


r.  (0) 


r. 

1 


Here  r.  ard  are  independent  variables,  but  r^(s)  is  a 
function  of  tinne  s.  These  fictions  trajectories  are  named 
trajectories  0.  Taking  t  so  that 

Ti  »  T  »T, 

we  make  the  average  of  each  term  of  Eq.  (2,  8)  along  the 
trajectories  0  over  s  from  0  to  t  .  On  consideration  of 

F^^^^(Pi,  ri(8);t)  F.^^^(Pj,  rj(s);t  +  s) 


m 


^  Fi^^^(Pi.r^(s);  t  +  8) 


m 


we  may  see  that 


t  ^  td-S 


(3.3) 


(3.4) 


(s);  t+  s) 
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=  J_  r^(—  +  —  .  1-1 )  r.(8);t+  8)  ds 

T  J  '  ;^8  m  ar.(8)  I  I  i 

A  1 


Pi  S 


+—  •  > 

'  at  m  3r^  i  '  0 


(3.5) 


where 


•Y 

<F.^^^X,;t)>  =  —  J  y  fJ^^  (p..  r.(8);  t  +  8)d8 
i  i  0  ^ 


By  coaree -graining  Eq.  (2,  8)  in  the  above  8en3e»  we  obtain 


(3.6) 


m 


i  7^ 


ar 


— )  <  F.^^^(X.;t)>  0  +  <J.>  »  0 

»  '  ‘  0 


(3.7) 


where 


J. 


1  y  ?i2  ._L  +^.  .-^  )  F..<^^X..X.;)dX. 

TT  ^  J  J  '  ij  ''  ^Pj  'J  '  J  J 


<^i)o  ' 


TV  j=ij^  ^  i  J  -J 

t  ■*  t+B 


4X.  ds 
J 


(3.8) 
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For  deriving  the  Boltzmann  collision  integral  fi*om  (3.  8), 
we  consider  the  following  three  points  which  are  commonly 
essential  in  statistical  investigation: 

1.  In  {3.  8).  the  summation  with  respect  to  j  is  made  first 
and  then  the  result  is  averaged  with  respect  to  s.  In  order  to 
make  the  time  average  separately  with  respect  to  each  of  particles 
j>  it  is  necessary  either  to  choose  t  sufficiently  short  that  particle 
i  interacts  with  particles  j  once,  at  most,  during  the  time  between 
t  and  t  -I-  T,  or  to  assume  that  the  interactions  are  very  weak. 
(Suppose  that  there  are  two  strong  collisions  of  particle  i  during 
the  period  between  t  and  t  T .  The  effect  of  one  collision  when 
it  occurs  first  is  different  from  the  effect  of  the  same  collision 
when  it  is  assume  d  to  occur  second.  Note  that  X^is  an  independent 
variable  and  is  the  same  in  the  two  cases,  but  in  the  first  case 

is  differaiU  from  in  the  second. )  The  necessity  of  the  above 

condition  will  be  discussed  again  with  respect  to  (3.  ZL). 

Z.  Since  we  are  observing  <  F.^^^(X.;t)>  ,  i  =  1,  2,  3, . 

1  i  0  ' 

it  is  impossible  and  also  unnecessary  for  us  to  know  the  precise 
time  when  the  precise  state  with  which  a  particular  particle  j 
begins  to  interact  significantly  with  particle  i.  Thus  we  may 
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set  the  same  assumption  as  the  Stosszahlansatz  by  Boltzmann 

with  respect  to  the  slate  with  which  particle  j  begins  to  interact 

significantly  with  particle  i  during  the  time  between  t  and  t  +  t; 

Particle  j  is  one  of  the  particles  which  are  located  near  particle 

(N) 

i  at  time  t,.  and  among  which  D  '  may  be  assumed  to  be  symmetric 

If  we  take  Ar,  for  the  domain  of  the  configuration  space  in  which 
D 

such  particles  j  exist ,  definition  2.11  leads  to; 


^(2) 

‘io 


tit+s 


ds 


(3.8) 


Note  that  the  domain  of  the  integration  with  respect  to  dry^  is 

limited  by  Ar^  in  accordance  with  the  definition  of  which 

(2) 

is  continuous  unlike  F,\  ;  Further,  on  coarse -graining  with 
respect  to  i  in  accordance  with  (2.  13),  we  have 

<~>0 


dX^  ds 

D 


t  t+8 


(3.8)” 


On  consideration  of  the  above,  Eq.  (3.  7)  yields 


<  +  ® 


F^^^(x^:t)> 


+  <  J  >  =0 


(3.7)' 
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3.  The  result  of  the  coarse>grainii;g  must  be  independent  of  t. 
Otherwise*  the  result  has  no  definite  physical  meaning.  The 
situation  depends  on  the  distribution  of  particles  in  A  r^  and  T. 

If  the  distribution  is  uniform,  the  result  is  perfect.  As  the 
density  of  a  gas  becomes  more  rar>ified  and  the  distribution  of 
particles  in  Ar^  and  t  becomes  more  non-nniform,  the  result  tends 
to  depend  on  T.  Then  the  possibility  of  kinetic  theory  diminishes. 

Instead  of  < J  >  ,  we  may  make  <J.'>  by  coarse-graining  J 
0  '  I  i 

over  time  from  t  to  t  +  t  along  trajectories  I  which  satisfy 


dr.(s) 

ds 


dr.(s) 

J 

ds 


?<(■) 


m 


II 

m 


dPi(s) 

ds 


-'^ii  <•> 


dp.(s)  ^ 

ds  J' 


(3.9) 


with  initial  conditions 


r.(0)  =  r.,  p.(0)  =  p. 


( Trajectories  I) 


r.(0)  3  r  .  p.(0)  =  p. 

J  J  J  J 


Thus  we  obtain 


r  T  \  I  C  T  f  ^  dp^Ce)  <) 

N-*!  >I  -  <-3r-  -J— tn  ‘ 


•^^j^(8).Xj(8);t+B)dXjd8 


(3.  10) 


(2.  11) 
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Equations  (3.  9)  constitute  the  characteristic  equations  of 
Eq.  (3.  1)  if  d/St  is  replaced  with  S/S  and  t  with  t  +  s. 

Therefore,  by  integrating  Eq.  (3.  1)  along  trajectories  I  given 
by  (3.  9),  we  obtain 

14-  '•if 

(  >J  T 

=  I  F..'  {x.(a).  X.(s);  t  +  s)!^  =  0  (3,  12) 


where 


d  _  S 

“  TT 


dr  (s) 


S 

S**.(8) 


(s)  9 

^  “"^i  ■  Sr .(s) 

J 


+ 


S 

dp-is) 


<ip.('^) 


s 

dpj(8) 


By  making  the  average  of  <  J.>  j  with  respect  to  particles  in  the 
same  way  s.s  is  shown  by  (2,  13),  we  obtain 

—  N  S  /Pb^*^  d  > 

J  J  ,  '"Ts - SCTs)  ds  •■3p^8)  ; 

o  ti  r.  a 


F^^^X^(s),  X^(s),  t  +  s)  dX^  ds 


(3.  14) 
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I 

I 


where  A  r.  is  of  the  same  implication  as  stated  with  respect  to 

D 

(3.  8)' .  On  definition  of 


r  d  ,  _  3  ,  Pa<‘>  3  ,  Pb'®’  3 

^  ds  de  m  *  dr^ts)  m  '  "Sr^Js) 


'dp  (s)  ^  ^ba^*^  dp.  (s) 


(3.  15) 


relation  (3.  14)  yields 


7  ’  -  N  f  f  f  r  **  1  .  ^ 

<J>  ■"TT"JoJArj^''-“Jir7^  Ts' 


Pa<«)  d 

m  *  dr^(8) 


'  d  12) 

•  JCT.)  >  '■■  '  *  'I 


_  (3.  16) 

Unlike  (3.  12)  [  ^ _ ]  F  in  the  above  integrand  does  not  vanish. 

/  ds  . 

ab 

The  reason  is  that 


Ssf^  -■ 


A  and  ^  considered  in  {Z,  13)  cover  domains  of  the 

configuration  space  which  are  larger  than  the  domain  of  the 

force  field  induced  by  a  single  particle;  hence  only  a  few  of  all 

the  particles  among  two  groups  of  particles,  one  group  in  states 

between  X  and  X  +  L  K  and  the  other  in  states  between  X, 
a  a  a  b 

and  Xj^  +  A  are  really  participating  in  binary  collisions 
of  the  mode  specified  by  (3.  13);  many  of  them  are  in  collisions 
of  different  modes  and/or  in  free  flight.  Hence  we  may  write 
for  (3.  16) 


<  J>  =  ^  t  ).  Xjj(T);  t  +  T  ) 


-TV  JoJir^  (JT  “■  ■^•1^ 


F*^*(X^(b),  X|^(b);  t  t  s)  dX  ds 

(3.  17) 

It  will  be  shov  that  the  first  integral  in  the  above  yields  the 
Boltzmann  collision  integral  while  the  second  integral  is  insignificant 


if  the  number  density  particle  is  as  low  as  usual. 

Particles  in  slate  may  interact  with  particles  in 
state  X  during  the  period  of  time  between  t  a^d  t  +  t  only  if 
r  at  time  t  is  ch'isen  in  the  domain  of  the  configuration  space 

D 

between  r'^  and  r^  respectively  defined  by 
(^■'b  ■  'a^  •  ^  Pb  "  Pa^  ^  ° 


f'b  •I’b  - 

I  Pi,  -  PaP"'  <*• 

as  is  illustrated  in  Fig.  3.  1.  Hence,  by  presenting  r^  with 
a  cylindrical  coorduiate  system  with  the  origin  at  r  and  the 

ol 

altitude  in  the  direction  of  p^  -  p^,  we  have 

i>  I  P^,“  P_l  <• 

^'b  =  J  ‘^'b  *  ■  - m - J 

(3.  19) 


Here  stands  for  the  direction  angle  and  h  the  radius  in  the 
plane  of  zero  altitude.  We  also  assume  that,  before  and  after 
a  collision  between  two  particles,  the  correlation  between  the  two 
particles  is  weak  so  that 


X^:t)  =  F'‘^{X^:t)  f'  MX^:t) 


F'^MX^{t).  X^{T):t  +  T  )  =  F'^^XjT);t+M  f'^^Xj^{t):  t  +  T) 


(3.20) 


On  consideration  of  (3.  19)i  the  second  member  in  the  right-hand 
side  of  (3.  17)  is  of  the  order  of 

N  f,r.»  ^  T\  i>  .-THr-T— - 


K  C  21^  i.  IIs  J 

■fV  J  '  «  -ir 


dp,T 


^  /VPal 


d  ^  h  dh 


K  d  r 


b  /  <‘  ^  V  ,,^1)  .  . 

-  ‘jr  •  -« — “  >  ^ 


(3.21) 


where  A  rj^  is  given  by  (3.  19)  and  is  the  volvime  of  the  configuration - 
space  domain  in  which  we  find  at  t  such  pa:  tides  as  lo  interact 
with  particles  of  state  X^  in  the  future  between  t  and  t  +  t.  It 
is  reasonable  ,  therefore,  to  assume  that 


0[  J  dtp  hdh],  =0  =5x10  cm'^ 


where  is  the  order  of  cross  section  or  o  is  the  maximum  distance 
between  two  particles  within  which  their  interaction  is  significant 


The  assumption  is  reasonable  only  for  neutral  molecules.  See 
for  example,  J.H.  Jeans,  The  Dynamical  Theory  of  Gases  (Dover, 
N.  Y.  ),  fourth  edition,  p.  ^81. 


and  that 

”1/3  ^  2 

(N/V)  «0[i  (p.  -  p,)/m]  «meam  free  path  length  =  — ^ - 

D  a  •>  /  _  _a 


2'  n  c=N/V 


By  choosing  t  as  in  the  above,  we  have 

^  NAr,  1 

i-rr—^  «  — rr —  «  — rz — 


=4=  0.  2 


(3.22) 


2/^  -2 

IfN/V  =  10  ,  then  o^(N/V)  '  '  =  5x10  ,  Under  this  condition, 

it  appears  barely  possible  to  choose  ^  so  that  (3*22)  is 
satisfied.  If  N/V  >  lO^'*,  then  =  5,  and 

(3.22)  does  not  hold  re-tardless  of  ^  .  It  is  noted  that 
condition  Vl^r^/V  I  is  the  very  condition  that  is 
short  so  that  particle  i  interacts  with  other  particles 
once  at  most  durln;  the  time  period  between  t  and  t  ♦  . 

Note  the  statements  which  follow  (5.8).  Under  condition  (3.22J, 
the  term  miven  by  (3.21)  is  nc-rli^ibly  small  as  compared  with 
the  first  two  terms  of  Eq.(3.7).  Thus,  with  the  help  of 
(3.19)  and  (3.20),  Eq.(3.17)  yields 


IPh  -  PJ 


d^  h  dh  d 


See  J.  H.  Jeans,  The  Dynamical  Theory  of  Gases 


. I . . 


are  of  microscopic 


If  we  consiier  thnt  T  .  r  (t')-r  and  r,  (T'^-r 

a  a  D  b 


scales  and  hence 


etc. 

in  the  time  scale  of  (  we  finally  have 

—  r. 


I  Pb-P, 


d  ^  h  dh  d 


(5.25) 


Since  *****  ^  "^i  V  I  made  by  averai^inij 

resoectively  alonw  trajectories  0  and  trajectories  I,  these 
two  averaf^es  are  not  mutually  equivalent.  In  the  followin;;, 
however,  it  will  be  shown  that 

after  their  beini?  coarse-grained  with  respect  to  :articles 
in  the  same  way  as  Is  suown  by  (2.13):  is  a  function 

of  and  t,  and  exists  along  a  line  named  in  the  X^t- 
space  as  is  shown  schematically  in  A  of  Fig. 5.2.  By  the 
assumption  oi  binary  collision,  line  coincides  with 
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one  of  th  -  trajectories  I.  Since  we  are  dealing?  with  a 
ainule  aystem,  vaniahea  unless  particle  i  is  in 
interaction  with  particle  j  so  tnat  J*.  .  is  significant. 
Almost  straight  lines  in  the  same  illustration  are 
trajectories  0»  In  B  of  Fig. 3*2,  ^“^i^  I 

are  shown  at  time  t.  anoears  as  presented  jy 

J -function;  *1®®  appears  as  if  J*-function, 

since  Jj  exists  only  on  one  of  the  tra je'.tories  I. 

On  the  other  hand,  made  by  averaging 

along  trajectories  0  is  shown  as  if  an  error  function. 
Obviously 


S^i>o  ^  <^i>i 

It  is  obvious,  however,  that 

r<ji(Xiit)X«i .  “i 

This  is  because  both  are  the  same  total  intensity  of 
from  t  to  t  +  '2"  divided  by  (  .  If  the  number  of 
particles  wuich  exist  in  states  between  and  X^  +  ^ 

is  sufficiently  large,  we  may  assume 


where,  of  course. 


*a‘  ^  *a 

<J.(X,,t)X  «i 

X 

a 


(5.29) 
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r  3  a 


<(ji(X,  ;t))i  dX, 


(5.30) 


Needless  to  say,  condition  (3»28)  is  compatible  with 
(J.26).  According  to  the  above,  Eq.(5.7)'  yields 


+ 


U 


U 


(3.31) 


In  order  to  reduce  Eq.(3*31)  to  Eq.(2.13)«  it  is 
necess  ary  to  replace  ^ 

(2.14)  and  (2.16); 


with  f  in  accordance  with 


f * (c,r;t> 


c  =  p/m 


(3.32) 

(3.33) 


It  may  be  a  trivial  matter  that  we  takv,  off  subscript  a, 
replace  subscript  b  with  1,  and  write  me*  and  mc^ 

respectively  for  P.j(^)  and  Pp('C^  completing  the 
reducti  on. 

Finally,  the  essential  points  in  the  aoove  derivation 
are  summarized  as  follows: 

1.  We  may  take  some  other  class  of  trajectories  in 
nlace  of  txajectories  0,  as  long  as  tney  are  equivalent 
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in  time  ec^le  t 

2.  Also,  trajectories  I  may  be  replaced  'vith 
trajectories  II,  etc.  ,  if  these  clas  es  are  equivalent 
in  time  scale  Tp. 

3.  ic  an  independent  variable,  and  Ff^^(X^;t) 

is  a  discontinuous  function  to  be  presented  by  S  -function. 

But  F^^^(X  ;t)  is  a  continuous  function,  if  N/V  and 
a 

^  X^  are  sufficiently  large.  The  situation  is 
similar  with  respect  to  ^  (X^,X^.;t ),  etc. 

4.  Althou'jh 

3),X.(s);tts)  =  0 

in  y«t 

fuf  Lb  f'  '>(Xa<b),X,(.),t*8)  f.  0 

in  (5.15). 

5.  The  domain  of  inteirration  ^  dr^^  considered 
in  (3.19)  is  finite  in  the  direction  of  p.-p  .  This 

D  3 

is  neres  -ary,  after  replacing  F^^^(X.  ,X  .;t )  with 

ij  '  J 

- 

F  (X  ,X,  ;t),  for  taking  into  account  the  fact  that 
a  b 

a  particle  with  its  finite  velocity  Is  localized  in  a 
narrow  domain  of  the  configuration  space  durln,'  a  short 
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period  of  time  Note  that  ,X.;t)  is 

^  J  ^  J 

dlscontlnuously  localized  by  itselfi  but  that 

F  ^^(X^,X^;t)  ie  a  continuous  function*  If 

were  symmetric  in  the  microscopic  sense,  the  above 
limitation  of  the  integration  with  respect  to  db^ 

would  not  be  .iustified*  Then  the  value  of  the 
collision  inte^rral  would  always  be  indefinitely  large. 


I 


’  I 

I 

i 


I 

I 

f 


t: 


IV.  DIFFICULTIES  IN  USUAL  DE'.^IVATIONS 


Thtre  have  been  differ^tnt  derivations  published 
by  seyeral  authors*  As  typical  examples,  the  derivations 
by  Kirkwood  and  Ross,  by  Bogoliubov,  and  by  Grad  are 
Investigated  in  the  following.  Because  of  the  assumption 
of  symmetry  of  considered  in  the  microscopic  sense, 

these  derivations  involve  common  difficulties. 


4.1.  Derivation  by  Kirkwood  and  Ross’^ 


1.  Being  coarse-grained  with  respict  to  time, 
Eq.(2*S)  yields 

Pi  n 


<,-f  * 


where 


.z 


F^^^X^;t^s)de 


(4.1) 


(4.2) 


By  assuming  that  the  distribution  is  symmetric,  we  may  replsce 


^J.  Kirkwood  and  J.  Rosa,  The  Statistical  Mechanical  Basis  of  the 
Boltzmann  Equation,  in  Proceedings  of  the  International  Symposium 
in  Statistical  Mechanics,  Brussels,  1965,  edited  by  I.  Prigogine 
(Interscience  Publishers,  Inc. ,  New  York,  1958. 


iO 


F^^^(Xi;t)  with  F'^^Xj^;t) 
Fjj\x.,Xj?t)  with  F^^^Xj^,X2;t) 


(4.3) 


with  N 


The  next  procedure  is  to  define  K' ‘  by 


F^^^(X^,X2;t+8)  = 


f 

K^^^Xj^,X2/X^,X^;8)  F^^^(X^,X'2;t)dX^  dX^ 


(4.4) 


Since  '  is  i^iveii  explicitly  by 


K^^^Xj^,X2/X|,X^;8)  =  J(X^-X£-4Xj^(s))  S(X^-^  -^X2(s)) 


(4.3) 


it  is  easily  shown  that 


F^^^Xj^,X2;t>a)  =  F^^\x3^-aXj^(s),X2  -iiX2(s);t) 


(4.6) 

According  to  Kirkwood  and  Ross.  ^X^(s),  standinf^  for4  P2^(s) 

and  A  rj^(8),  and  AX2(8).  standing  for  4  P2(s)  and  4  r2^®^’ 
are  the  changes  in  momenta  and  spoce  coordinates  of  particle  1 
and  particle  2  determined  by  solving  the  equations  of  motion 
of  the  two  particles  in  accordance  with  the  as  umption  of 


binary  collision: 


d  A  PjCo)  ♦  APj^Ca) 


d  s 


m 


and 


dAp^Cs) 

-r:--  = 


(4.7) 


similar  squaticns  for  particle  2  . 

Under  this  circumstance,  they  claim  that  given  oj  (4.6) 

is  a  solution  of  Eq.(3.1)  where  ^/c)t  is  replaced  mithc^/^e: 

(ji  ♦  Ji.i.  .  Ji  .  ±  ,  .  -i.  ♦  0 

^s  m  jj^rj^  m  ^r2  12  ^  Pj^  ^  ^^2 

(4.8) 

Here  s,  p^,  p.*  r^  and  r^  are  to  be  independent  variables. 

It  is  obvious,  however,  that  given  by  (4.6)  under 

conation  (4.7)  does  not  satisfy  £q.(4.8).  Instead, 
it  is  obvious  that  satisfies  Eq.(4.8),  if  a!  X. (s) 
and  X^{b)  are  solutions  of 


d4ir^(s)  Pj^  ^AP^(8) 


d  s 


m 


ds 


-J 


12 


(4.9) 


etc. 


in  place  of  (4.7).  Note  that  ?2*  *’i  **2 

Indepenc'ent  of  s  in  (4.9),  and  h>?nce  t;.e  equations  are 


J2 


not  the  equations  of  notion  of  the  two  particles*  Only 
when  ^  X^(s)  and  ^X^Cs)  are  extremely  small,  li:qs*(4*9.) 
coincide  with  Eqs.^4.7)*  It  is  noted  that  the  particular 
solution  of  Eq.(4.8)  Integrated  along  trajectories  determined 
by  Eqs*(4.7)  is  considered  in  (3*12)* 

2*  Although  ^X^(s)  and  a  X.^sb)  satisfying  Eqs*(4.9} 

are  not  the  proper  trajectories  of  the  two  particlds  in 

interaction,  such  as  expected  to  be  by  the  authors,  yet 
(2) 

F  given  by  (^.6)  is  a  correct  solution  of  Eq*(4.8) 

under  condition  (4*9)*  Therefore,  if  the  succeeding  treatncmt 

is  proper,  it  may  lead  to  a  correct  result.  By  substituting 

(?) 

P'  '  given  by  (4,6)  in  the  collision  integral  of  (4*1), 
however,  we  will  meet  a  great  difficulty  in  the  following 
trcatrjnt:  Since  X^  and  X^  are  independent  variables,  the 
domain  of  the  integration  with  respect  to  dr^  is  to  be  the 
entire  configuration  sp^kce;  snd/or  ^  are  functions 

of  and  ,  and  are  completely  independent  of  ^r^(s)  and 

A  r^(s)  by  (4*9),  Therefore,  is  not  a 

factor  which  affects  the  magnitude  of  ^^2*  when 

extremely  large,  yet 

vanish  if  r^  -  r^  is  small,  and  vice  versa.  Hence, 
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tnere  ir  no  reason  for  limiting  tne  domain  oi  inte  r.iwion 
with  respect  to  dr_  in  the  direction  of  p^  -  o,  by 
■^jP2-Pi|/®  that  is  a  function  of  . 

3*  The  d\ f ’'iculty  mentl oned  in  tne  above  seems  to 

have  stemmed  from  the  authors'  belief  that  the  X's  and 

t  must  be  always  independent  variables.  This  belief  is 

(H) 

a  rational  consequence  of  their  assumption  that  J 

ven  in  (2.11  is  always  symmetric  with  respect  to  tne 

interchange  of  coordinates  between  two  similar  particles. 

Indeed,  accordant;  to  the  assumption,  it  is  pointless  to 

say  that  an  individual  particle  moves  alonr;  a  continuous 

trajectory;  there  is  no  such  individual  particle  as 

distinguishable  from  the  others  by  its  continuous  trajectory. 

In  spite  of  tne  technical  difficulty  mentioned  thus  ''ar, 

it  appears  that  the  authors'  approach  is  at  least  persistent 

(ID 

accordin'?  to  their  belief  th«t  D  is  symmetric,  but  tne 
situation  is  not  necessary  sc;  ^esinst  the  aut.iors' 
expectation,  ^Xj^(o),  X^  -^X^Cs);!)  ,;iven  by  (4.6) 

is  no  more  symmetric  as  is  snown  below.  Suppose  that  we  corn  ier 
at  y,  =4X,(s)  and  X^  X-(s).  Tnen,  accordin'  to 
the  symmetry  assumption,  it  must  nold  tnat 

F^^^(0,0;t)  =  X^(s),  d  Xj^(B)-4X2(«);t) 

Since  A  X^(s)  -  ^X^(s)  may  increase  as  s  increases, 
the  above  eouation  implies  that  (X^  ,X2;  t )  is  uniform 

alone  a  line  given  by  Xj^  =  If  disol.ice  the  origin 

of  the  coordinate  system  by  keeping  the  axes  as  parallel  to 
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those  before  the  di splacementf  we  may  realize  tnat 

(2) 

F  is  uniform  (Invariant)  ever  the  enitire  phase  space. 
This  obvious  paradox  is  caused  by  the  attempt  to  represent 
a  microscopic  collision  in  terms  of  which  is  assumed 

to  be  symmetric. 


4.  Microscopic  collisions  are  conceivable  only  when 
we  bet;in  a  theory  by  taking  D'  of  a  single  system  as  given  by 
(2.3)  which  is  not  symmetric.  By  taking  =  X^(s), 
we  are  observing  D'  along  the  trajectory  of  particle  i. 

It  is  worthwhile  to  note  that  there  are  various  solutions 
of  the  partial  differential  equation  (4.8),  according 
to  our  choice  of  functions  of  s  to  be  assigned  to  p^.  p^t 
rj^,  and  r2; 


Pt  »  pj(8),  =  rj(s) 

P2  =  P^(s),  r2  =  r*(s) 


(4.10) 


In  the  special  case  where  these  functions  given  in  the 
above  satisfy  the  characteristic  equatiens  of  £q.(4.8) 
which  are  the  very  equations  of  motion  given  by  (3*9) « 

is  invariant  with  respect  to  s.  Otherwise, 
is  not  invariant. 
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In  view  of  the  above  investigation,  it  appears  that  Kirkwood  and 
Ross  believed  the  symmetry  of  in  the  microscopic  sense. 

(Another  evidence  of  their  belief  is  found  in  ore  of  Kirkwood' s 
previous  papers  also.  )  This  belief,  in  spite  of  their  reasonable 
approach  of  coarse-graining  with  respect  to  time,  seems  to  have 
jeopardized  their  derivation.  In  general,  the  same  difficulty  as 
is  found  in  the  derivation  b>  Kirkwood  and  Ross  is  pointed  out 
commo.-'ly  in  derivations  beginning  with  the  symmetry  assumption 
of  o'”'. 

r 

4.  2.  Derivation  by  Bogoliubov 

The  derivation  by  Bogoliubov  is  similar  to  the  derivation  by 

Kirkwood  and  Ross.  Bogoliubov  defined  "streaming  operator"  which 

12) 

is  of  the  same  implication  as  that  of  K  given  by  (4,  5).  Therefore 
the  difficulty  seen  in  Bogoliubov' s  derivation  is  of  the  same  nature  as 
that  of  Kirkwood  and  Ross. 

♦ 

4.  3.  Derivation  by  Grad 

The  derivation  by  Grad  is  anumolcus  in  the  sense  that  the 

(N) 

derivation  is  consistent  in  spite  of  the  symmetry  of  D  assumed  in 
the  microscopic  sense.  This  paradoxical  situation  is  explained  in  the 
following.  Grad  defined  "truncated  distribution  functions",  e.g.  ,  Fj^(Xj) 


4 

J.  Kirkwood,  J.  Chem.  Phys.  14,  180(1946).  Particularly  the  paragraph 
beginning  in  pg.  181  and  ending  in  the  next  page  is  of  interest  from  the 
present  viewpoint.  For  a  system  to  be  investigated,  he  chooses  one  at  undom 
among  those  systems  constituting  a  canonical  ensemble.  As  a  result,  D'  '  of 
the  system  is  statistical  from  th  '  beginning  of  his  theory. 

^N.  N.  Bogoliubov,  Problems  of  a  Dynamical  Theory  in  Statistical  Physics 
(1946),  translated  byjiJ.  G.  uora,  in  Studies  in  Siaiisucal  Kaechamcs,  edited 
by  J.  de  Boer  and  G.  E.  Uhlenbeck(North-Holland  Publishing  Co.  1962)  Vol.I. 

$ 

H.  Grad,  Principles  of  the  Kinetic  Theory  of  Gases,  in  Handbuch  der  Physik 
edited  by  S.  Flugge  (Springe r -V erlag  ,  Berlin  1958),  Vol.  XII  pp.  205-295. 


36 


I 


I 


t 

s 


is  the  expectation  of  finding  no  molecule  within 
a  distance  of  molecule  1«  granting  that  molecule 

1  ia  in  state  X^.  It  is  also  assumed  that  the  sphere  of 
radius  (J*  limits  the  domain  of  influence  of  molecule  1. 

Qrad  derives  the  equation  ox  evolution  of  F^(X^)  from 
the  Liouvllle  equation*  instead  of  the  equation  of 
evolution  of  This  is  equivalent  to  considering 

the  distribution  function  ox  molecules  as  in  states  of 

free  flight  in  the  same  way  as  is  intended  by  the  Boltzmann 

a 

equation  .  Also  by  assuming  the  existence  of  the  sphere 
of  Influence,  the  further  treatment  appears  as  if  a  treatment 
of  collisions  among  elastic  spherical  bodies 'of  diameter  (T', 

instead  of  collisions  by  means  of  force  fields.  Therefore 

% 

the  spatial  domain  ot  the  collision  Integral  with  respect 
to  a  field  particle  is  always  the  space  cutside  the  sphere 
of  influence  of  particle  1.  In  this  way,  Qrad  could 
avoid  the  cause  of  the  common  difficulty,  that  Is^  the 
conslder/ttion  of  trajectories  of  particles  during  their 
collisions.  In  order  to  complete  the  derivation  by 
giving  the  equation  obtained  above  the  same  implication 
as  that  of  the  Boltzmann  equation,  however,  it  is  necessary 
to  assume  the  existence  of  "molecular  order"  in  the  behavior 


*  0  T 

Fj  (Xj)  is  equivalent  to 

t+s)  ds 

if  collisions  are  rare.  ° 
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of  two  particles  within  the  distance  O'  >  This 
assumption  is  not  compatible  with  the  postul  it i on 

/u  \ 

that  the  symmetry  of  v  and  hence  the  Bogoliubov- 

Born-Qreen-Kirkword-Tvon  hierarchy  are  valid  in 

tl.e  microsco'*lc  sense.  This  condition  is  out  of 

the  sc:.pe  Oi  Qrad's  theory.  Grad's  scheme  is 

ir.geneous  in  the  sense  that  it  has  separated  itself 

from  the  treatment  of  collision  mecnanism  wuich  is 

(V) 

not  compatible  with  che  symmetry  assumption  of  D 


V.  CONCLUSION 

1.  The  feasibility  of  the  Boltzmann  equation 
for  neutral  molecules  has  been  tested  and  confirmed 
by  experiments. 

2.  The  Boltzmann  collision  integral  is  derivable 
from  the  Liouville  equation  only  by  considering  that 

(n) 

the  distribution  function  D'  of  a  system  is  not 
symmetric  in  the  microscopic  sense. 

3.  The  microscopic  symmetry  of  P'  \  by  which 

a  particle  is  indistinguishable  from  the  others,  is  not 
compatible  with  the  Boltzmann  collision  integral  where 
a  particle  is  distinguishable  during  its  collision  by 
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its  continuous  and  ordsred  trajectory* 


4,  By  accepting  the  Boltzmann  equation  as 
feasiblej  we  must  admit  that  kinetic  theory  is  to 
begin  with  ^  which  is  not  symmetric* 
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F'.g.  2.2.  In  the  Boltzmann  equation,  the  effect  of  a 
collision  ie  an  integral  over  the  entire  phase  of  the  collision  as  is  shown 
in  A.  On  the  other  hand,  the  effect  of  a  collision  in  Eq.  (2.8)  is  accounted 
for  only  at  X.  between  X  and  X+dX  as  is  shown  in  B.  Reduction  of  the 
Boltzmann  equation  from  Eq.  (2.  8)  is  made  possible  by  coarse-graining 
Eq.  (2.  8)  with  respect  to  particles  and  time. 
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1  J  A  H  *  C  T 

The  derivation  of  thf  Boltzmann  equation  from  the  L'.ouville  equation 
can  be  a  means  of  testing  the  feasibility  of  a  kinetic -theoretical  '.'cheme,  since 
the  Boltzmann  equation  is  known  as  empirically  plausible.  From  this  view¬ 
point,  we  investigate  kinetic -theoretical  schemes  deriving  the  Boltzmann  jquation, 
We  may  derive  the  Boltzmann  equation  consistently  by  considering  the  Liouville 
equa  ion  to  govern  the  distribution  of  a  system  in  the  phase  space  which  is  not 
symmetrica  in  tl.e  microscopic  sanse.  On  the  other  ^nd,  we  cannot  derive  the 
equation  consistently  on  the  assumption  of  symmetry  of  the  distribution;  the 
Boltzmann  collision  integral  is  not  <_ompatible  with  the  symmetry  assumption. 
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